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Abstract 

The fermionic dispersion relation in the presence of a background magnetic field and 
a high temperature QED plasma is calculated exactly in the external field, using the 
Hard Thermal Loop effective action. As the field strength increases there is a smooth 
transition from the weak- field (qB <C q 2 T 2 ) thermal dispersion relations to the vacuum 
Landau levels when the background field is much stronger than any thermal effects 
(qB 3> q 2 T 2 ). The self-energy at finite field strength acquires an imaginary part. The 
spectral width becomes important for critical field strengths (qB ~ q 2 T 2 ), necessitating 
the use of the full spectral function. It is shown that the spectral function satisfies 
the usual condition of normalization and causality. Using the exact spectral function I 
also show that the production of chirality in an external electromagnetic field at high 
temperature is unaffected by the presence of the thermal masses of the fermions. 
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1 Introduction 

During the last few years an increasing understanding of high temperature gauge theories has 
been achieved through the consistent resummation of leading diagrams called Hard Thermal 
Loops (HTL) H|. The HTL resummed theory has been used to improve IR divergences in 
a number of processes, the most famous one being thermal gluon scattering and the related 
gluon damping problem ||. The HTL effective action has been constructed in several different 
ways, and it yields a gauge-invariant extension of the electric screening mass for non-static 
fields. Since strong background fields are likely to have been present in the early Universe it is 
of interest to look for solutions to the effective HTL equations of motion in such backgrounds. 
Most of the literature so far has been concentrated on the small fluctuations around the trivial 
background with zero field. It has been known for a long time that static magnetic fields are 
not screened by the HTL resummation and thus a constant magnetic field is also a solution to 
the resummed effective equation of motion. In this paper I explore the fermionic excitations 
around such a constant field at very high temperature. One interesting point is that this can 
be done exactly, i.e. to all orders in the external field. At the end we recover two well-known 
limits. On the one hand there is the zero-field limit where the thermal dispersion relations 
are well known 0, [5|. On the other hand, when the field is very strong compared with the 
temperature, thermal corrections become less important and in that limit we recover the 
zero-temperature Landau levels. For intermediate field strengths the self-energy acquires a 
non-negligible imaginary part due to synchrotron radiation and scattering with the heat bath, 
even above the light-cone. It is, therefore, necessary to study the full spectral function and 
not only an on-shell relation for the real part. It turns out that in the lowest Landau level 
the spectral width is rather large, when the field strength gets comparable with the thermal 
mass squared (qB ~ .M 2 ), and the quasi-particle picture is not reliable. 

Since there is a mass gap in the thermal fermionic spectrum even for a chirally invariant 
theory, it is not immediately obvious how the chiral anomaly mechanism works. The standard 
level crossing picture is not applicable since no levels ever cross the Dirac surface. I show, using 
the full spectral function in a background of electric and magnetic fields, how the spectral 
weight associated with particles and antiparticles can move continuously between the positive 
and negative energy solution without crossing the Dirac surface, and through this mechanism 
satisfy the anomaly equation. 
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In section |2|, I describe the method of diagonalizing the HTL effective action in a back- 
ground magnetic field, and compare the exact result with an approximate weak field formula 
in section |3|. The spectral function is calculated in section |4j. The main issue of the paper, 
namely the anomaly mechanism at finite temperature, is treated in section |5| for the case of a 
free field in 1+1 dimension and in the HTL approximation in 3+1 dimension. Some properties 
such as normalization and causality are discussed in an appendix. 

2 Dispersion relations from Hard Thermal Loops 

The HTL effective action for QED can be written as |J 

+ ^(H-m^-A^^— , (2.1) 
where = id^ — qA^ and the average (•) is defined by 

(/(«o,tt)> = , (2.2) 

where uq = 1 and u is a spatial unit vector. The thermal mass of the photon A4 2 is given 
by q 2 T 2 /9 and for the electron we have M. 2 = q 2 T 2 /8. The equation of motion for \1/ that 

is 

r / \ i 

* = . (2.3) 



follows from Eq. (|2.1 



/ u 

Equation ( |2.3| ) is a non-local and non-linear differential equation, which is, in general, very 
difficult to solve. What makes this equation much less tractable than the thermal Dirac 
equation, in the absence of an external electromagnetic field, is that the average over u is 
difficult to perform explicitly since [ILj, n„] = —ieF^ u ^ 0, i.e. not all components of n M can 
be diagonalized simultaneously. We shall in this section only deal with an external magnetic 
field and fix it to be in the ^-direction. The solutions to Eq. ( |2.3| ) in vacuum (Ai e = 0) 
are given by the standard Landau levels. Since the spatial symmetries of the system are 
unchanged by the thermal heat bath, we expect the eigenfunctions to have the same spatial 
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form as at zero temperature. In fact, after performing the u- integral in Eq. ( |2.3| ) the result can 
only be a function of the invariants IIj_, pi and p% and the 7-structure has to be proportional 
to 7IIL, 70P0 and j z p z $\ 

We shall therefore compute the matrix elements 



it'- 



ll ■ n 



(2.4) 



between the vacuum eigenstates. To be specific we use the gauge = (0, 0, — Bx, 0). Then 
the eigenstates are given by 



= exp[i(-pat+p y y+p z z)]I n . iPy (x) 



(2.5) 



In;p y {X) 



■\qB\' 



1/4 



7T 



exp 



v n\ 



2\qB\ [x 



Py_ 

qB / 



where k = {po,n,p y ,p z } and H n [x) are Hermite polynomials defined by 

ff„N = (-ireTA_e-V . 



(2-6) 



(2.7) 



These states form a complete set of functions in four dimensions when the energy is off shell. 

In the chiral representation and with qB > 0, suitable spinors can be formed from $ K as 

*f? K = diag[$ K , $ K , $ k _i]xk where Xn is a space-time-independent spinor, which can 

be determined from the Dirac equation. The vacuum Dirac operator in Eq. (|2.3|) gives by 

construction an eigenvalue when acting on \I/ K , but it is more difficult to determine the action 

of the thermal part since $ K cannot be an eigenfunction to u-U for all u. One way of calculating 

the matrix element in Eq. (|2.4j ) is to find a basis such that v ■ H\v p ) = v ■ p\v p ) and insert a 

unit operator J[d A p}\v p )(v p \ into Eq. (TO) ). (We use the notation J[d n p] = J d n p/(2ir) n .) The 

1 We use the notation a ■ b± — a x b x + a y b y for any two four-vectors a and b. In our convention three- 
vectors such as p — (p x ,p y ,p z ) and 7 = {^ Xl ly,lz) are the contravariant components of the corresponding 
four- vector and thus have Lorentz indices i = 1,2, 3 upstairs, i.e. p x = p 1 etc. We use the Minkowski metric 
diag(+, — , — , — ) so that p% = —p l and 7; = —7* for i = 1,2, 3. 
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unit operator is, of course, independent of v, so in particular we can choose v = u and change 
the order of integrations between p and u. In the gauge we use, an eigenvector to v • II is 
given by 



(x\v p ) = exp 



-ip t + ip z z + ipyy + i p x x + 



qBvy x 2 * 



[2.8) 



After computing the matrix elements in Eq. (|2.4|) we find indeed that they are diagonal in k 
for u and u z , and have a mixing with the first subdiagonals for u x and u y . We define (wo,2,±) 
by 



(2-9) 



(2.10) 



and k =F 1 = {po ; n T ^-iPyiPz}- These are exactly the components that occur naturally when 
we include the 7-matrices in the chiral representation. The explicit calculation of (wo,z,±) is a 
bit lengthy but straightforward and is done by performing the integrals over x, x', p and u in 

{<P Kl \x'){x'\u v } — 

p 



($ K ,\(^-)\$ K ) = fdxdx' dp^($ K ,\x)(x'\u p )^(u p \x)(x\<S>J . (2.11) 
\u ■ 11/ J Air u ■ v 



The result reads 



(u )k 



1 



n\ V27T ^-00 

Pz , Po+P 



X 



In: 



dsH 2 n (s)e- s2/2 

p Sy / 2qB 



Sy/2q~B 

2p 2 ~ Po-p z ' 2p 2 Jp\ - p 2 ~ ~~ 2Jp 2 - p : 



+ 



arctan • 



:2.i2i 



i= r dsH 2 n {s)e-° 2 l 2 

/ Z7T J —00 



n\V2n 



Po(2p 2 z - qBs 2 ) Po+p 



4 p 4 



In 



p z (2p 2 -p 2 ) s^2qB 
H — ; , arctan 



Po-Pz 
Sy/2q~B 



2p 4 



Po~P 2 



2 VPo - r 



(2.13) 
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(«+)« 



dsHJs)H n -i(s)e 



-s 2 /2 



X 



2im\{n - 1) 
Sy/2qB p Q sp z y / 2qB Po+p z 



\ 2p 2 2p 4 
2p 2 z (p 2 -p 2 )-p 2 qBs 2 



In 



arctan 



s^2q~B 



<«->« 



;2.14) 



;2.i5) 



where p 2 = p 2 + qBs 2 /2. With these definitions the Dirac equation effectively reduces to a 
4x4 matrix in the spinor indices, since the other quantum numbers have been diagonalized. 
In the massless limit (m = 0) the left- and right-handed parts factorize and the Dirac equation 
takes the form 



fl-M 2 elfl 



u ■ n 



x 





D r (k) 



D l {k) 








(2.16) 



where 
D r (k) 

dm 



-Po + Pz + M 2 e ((u ) K - (u g ) K ) 
-t^BE-M 2 e (u + ) K 

-Po ~Pz + M 2 e ((u ) K + (u z ) K ) 
i^/2q~Bn + M 2 e (u + ) K 



i^/2qBn - M 2 e {u_) 
-Po ~Pz + M 2 e ((u ) K -i + 



-W2qBn + M 2 e (u-) K -i 

~P0 +Pz+ MI{{uq) k -X - (u z ) K -i) 



(2.17) 



In the lowest Landau level (n = 0) Eq. ( j2.16j ) reduces to a 2 x 2 matrix, since only one 
orientation of the magnetic moment is possible. It is easy to take the determinant of Eq. ( |2.17| ) 
to find the dispersion relations, which for the right-handed component are 



n > 1 : 



n = : 



(po ~Pz- M 2 e ((u ) K - (u z ) K )J 

X (po +Pz~ MI{{Uq) k ^ + (U Z ) K -!)) 

-y2q~Bli-iM 2 e (u + )^ 2 
Po-Pz- M 2 e ((u ) K - (u z ) K ) 



X18) 
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Figure 1: Dispersion relations for the right-handed branch in the lowest Landau level 
(n = ), neglecting the imaginary part. As the B-Eeld increases thermal effects become 
less important and the dispersion relation approaches the light cone, which is indicated 
by solid lines. All dimensionful parameters are given in units of the thermal mass A4 e . 



These relations are only meaningful for stable propagating quasi-particles with well-defined 
relations between momentum and energy, i.e. when the imaginary parts are negligible. In 
general there are imaginary parts in the functions {uq jZi ±) k , which are discussed in Section [|. 
It is anyway useful to first solve Eq. ( |2.18| ), ignoring for the moment the imaginary part, since 
the zeros of the real part indicate where the spectral functions are peaked, at least when the 
imaginary part is small enough. This can conveniently be done numerically as all the integrals 
in Eqs. fl2.12|) to (|2.15 ) are well convergent. The dispersion relations for several field strengths 
in the lowest Landau level are shown in Fig. ^ 

In the lowest Landau level a right-handed particle (positive chirality, x — R — +1) with 
q and B positive, can only propagate in the positive z-direction since the magnetic moment, 
and thus the spin, has to point along the field. Positive chirality implies positive helicity for 
particles and thus positive p z . The left-hand side of Fig. [I] has to violate one of these sign 
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arguments. The solution is that hole excitations have opposite chirality-helicity relation and 
can thus propagate for negative p z . As the field increases the hole branch develops a new 
sub-branch and disappears continuously for large enough fields. The new branch must not be 
taken too seriously, since it only appears when the imaginary part is non- negligible and then 
only the full spectral function is meaningful. The particle branch approaches the light cone, 
i.e. the vacuum dispersion relation, in a smooth way as the field strength increases. This 
is physically very reasonable since, for very strong field strengths, the thermal effects should 
disappear. Once again it should be emphasized that the above analysis is based only on the 
real part of the self-energy and it can only serve as a guiding line to describe what kind of 
modes propagate in the plasma. For a complete description, which is necessary for qB ~ M.\, 
where the imaginary part is comparable with the real part, the full spectral function has to 
be used, as we do in section f§. 

The HTL effective action is derived under the condition that the temperature is much 
larger then the momentum. Here, the magnetic field enters only through the covariant mo- 
mentum and should thus satisfy the condition II 2 ~ qB -C T 2 . On the other hand, already 
when qB ^> M. 2 e ~ q 2 T 2 (which is the only scale where T enters in the HTL approximation) 
the thermal corrections start to get small compared with the tree- level part. Thus, for small 
coupling the HTL corrections become small before they are invalid. 



3 Comparison with an approximate formula 

In a direct one-loop calculation of the fermionic self-energy in a magnetic field |7j one is 
naturally led to an approximation where the full 5-dependence is kept only where it is added 
linearly to the momentum squared. In other places it enters only to 0(B 2 ) (see |7j for details). 
Since the result from this approximation is surprisingly simple, it is worth commenting on its 
relation with the exact solution. From |7| we find the Dirac equation 

[ji- m -E(p 0s p^n ± )]* = 

s(p Q ,n 2 ) loPo -r(p ,n 2 ) lzPz -r( Po ,(n-j) 2 )fl ± -m\ = , (3.1) 
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where II 2 = II^+p 2 . The functions s(po, II 2 ) and r(po, II 2 ) are derived from the HTL effective 
action without background field and they are given by: 



p s{p , n 2 



p z r(p ,U 2 



Po-M 



2 / U 



u ■ p 



Pz- M 



u ■ p 



M 2 
Po r-^r hi 

^° 2 n 



Po + 


|n| 


Po - 


|n| 



(3.2) 



Pz 



n 2 I 



Po 

2 in 



In 



Po + 


\U\ 


Po - 





(3.3) 

It is almost possible to guess the expression in Eq. ( |3.1| ) from the standard expression for 
the HTL Dirac equation [|, |]. The usual momentum p^ should be replaced with the gauge- 
invariant momentum II M , but there is an ambiguity in replacing p 2 by II 2 or by HH. The 
correct way follows from the calculations in [|7|]. 



The difference between Eq. (|3.1| ) and the exact formula is related to the order of doing 
the average over u and replacing p by II. Comparing the exact expression for (uq) with 
Eq. ( |3.2|) , before doing the ^-integration, we would like to specify under which circumstances 
the approximate equality 

,dtt 



nW2n 



dsH 2 (s)e- s2/2 - 



4ir 



dn 



u 



47T u po - u z p z -u-p± 



up 

UqPq - u z p z - 



(3.4) 



is valid. The first term comes from the exact expression and the second from the approximate 
formula Eq. ( |3.2| ). In a standard coordinate system with u = (sin 9 sin (f), sin 9 cos 0, cos 8), so 
that u± = sin9, we see that the difference lies in the integral over the azimuthal angle (f). 
In the exact formula the integral over is replaced by a more complicated integral over s, 
involving the exact external states, i.e. the Landau levels. The two expressions in Eq. ( [3.4|) 
do not coincide, except in some particular limits. Expanding both sides of Eq. ( |3.4|) formally 
in powers of qB and LT^, we are led to comparing the integrals 

^B\ 2k ..2„, mfc (2fc)! min( ^ 



1 



2,71 J— 



ds H n (s)e 



-s 2 /2 




(uiqBY 



2 2k 



E 

z=o 



2 l n\ 



(n - l)\(k - 
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- + , (3.5) 



77" («± ■ P±J 
Z7T 



where we used the fact that = qB(2n + 1) when acting on a Landau level. First we notice 
that for k — 1 the two integrals coincide. Then, we find that the leading terms in the limit 
of large n, for fixed k, also coincide. We can thus expect that the approximative formula 
Eq. fl3.1|) is useful both for weak fields and for very high Landau levels. It should, however, 
be noticed that the expansion converges badly close to the light cone, and that it eventually 
breaks down for hole excitation of high momentum JF]. 

In many applications it is only the dispersion relation for small momenta that is important. 
Using Eq. (|3.1|) in the limit m = 0, we can easily obtain an approximate formula for the 
dispersion relation in the lowest Landau level around p z = in the presence of a weak 
magnetic field and x — 1 

w^+SdM'-m) ' (3 ' 6) 

For the right-handed branch of the lowest Landau level in a weak magnetic field (qB = 
0.2 .M^), the dispersion relation following from Eq. ( |3.1[) is shown in Fig. where it is also 
compared with the exact solution of Eq. Q2.18| ). 



4 Spectral function 

The dispersion relation was solved in section |2], ignoring the imaginary parts of the self-energy. 
This is only a good approximation for small magnetic fields, where the imaginary parts are 
small. Since we have the exact expression for the self-energy it is not too difficult to study 
the complete spectral function and to see how good the quasi-particle picture is. The spectral 
function can be defined from the representation of the retarded and advanced propagator as 

s(£±i£ , p) = £^^§|L . (4.i) 

In the real-time formalism of thermal field theory the spectral representation goes through in 
much the same way as at zero temperature, the only essential difference being doubling the 
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Figure 2: Comparison of the dispersion relation from the HTL effective action and the 
weak held approximation in the lowest Landau level for qB = 0.2_A4g. All dimensionful 
parameters are given in units of the thermal mass A4 e . 



degrees of freedom (for a recent review see H). The full thermal propagator takes the form 
i 

S^(E,p) = r dE'A(E',p)cr z B-\E') ( ) B(E') , (4.2) 

J -°° \ U E-E'-ie ) 

where cr z is a Pauli spin matrix and 13(E) can be chosen to be 

13(E) =( {e [e x + . (4.3) 



For a free Dirac fermion we have 

A(E, p) = ('JqE - 7P + m)sign(E)5(E 2 - p 2 - m 2 ) , 



(4.4) 
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but in general A(E, p) can have both a 5-function part for the quasi-particles and a continuous 
part. The HTL fermion propagator without any external S-field is given by 



S(E,p) 



s(E,p)^ E - r(E,p)jp + m 
s(E,p) 2 E 2 — r(E,p) 2 p 2 — m 2 



(4.5) 



where the functions s(E,p) and r(E,p) are defined in Eqs. ( |3.2| ) and ( |3.3| ). For E > \p\ the 
only imaginary part comes from the analytic continuation using ±ze and the contribution to 
the spectral function become 5-functions at the solutions of the dispersion relation. Below the 
light cone, i.e. when E < \p\, there is a finite imaginary part emerging from the logarithms in 
Eqs. (|3.2j ) and (|3.3| ) giving a continuous contribution to A(E, p) . In the appendix it is shown 
that S(E,p) fulfils the general requirements of normalization and causality. 



4.1 Spectral function in the presence of a 5-field 

Since the self-energy in the presence of the 5-field does not have any singular points away from 
the real axis, and since the HTL corrections are negligible for large complex E, we expect 
that the propagator still has the correct analyticity properties and that the normalization 
and causality properties, discussed in appendix [A|, are satisfied. We have checked the sum 
rule in Eq. ( |A.3| ) by direct numerical calculations and it is indeed satisfied. The analytic 
continuation E —>■ E ±ie is more complicated in the presence of the background field, but it 
can be summarized by the formula: 



1 



: arctan 



s^[2qB 
2^W 



9(E 2 -p 2 



arctan 



^E 2 - p 2 
+ 6(p 2 -E 2 )i [ e(p 2 z -E 2 ) 

+ 6(E 2 -p\ 



2^W 
1 



2^p 2 - E 2 
sign(s x 



p- 



ln 



2^p 2 — E 2 — s^2q~B 



2^p 2 — E 2 + s^2q~B 



In 



2^p 2 - E 2 

i7rsign(s)sign(_E) 
2^p 2 - E 2 



\s\ y/2qB 


-2^/p 2 


-E 2 


\s\ yJ2qB 


+ 2^ 


-E 2 



(4.6) 
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Figure 3: Spectral function for various magnetic field strengths at the momentum 
p z = 0.0 in the lowest Landau level. For weak fields there are two equal peaks around 
E ~ ±1.0 (indicated by the vertical solid lines), corresponding to the thermal mass 
of the particle and anti-hole solutions. As the field increases the width of the peaks 
increases and the positions are shifted to a slightly higher value. For intermediate fields 
the spectral function is very wide and eventually it gets more concentrated around E = 
0.0, which is the position it should have without thermal correction. All dimensionful 
parameters are given in units of the thermal mass M. e . 



for E -> E - ie. 

With this explicit expression it is straightforward to calculate the spectral function nu- 
merically. We have done so for p z = and p z = 0.5 for the right-handed branch in the lowest 
Landau level; the result for the term proportional to 70, namely 



A? LL (E, Pz 



tr 



-(1 + 75)70^, Pz,n = 0) 



= ^-(S? LL (E-ie,p z )-SK LL (E + i6, P ; 



S£ LL (E,p z 



p -p z - M 2 e ((u ) - (u z )) 



(4.7) 
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is presented in Figs. [| and |j. These figures should be compared with the solution of the 
dispersion relation from the real part in section and correspond to two vertical cuts in 
Fig. |I] at p z = and p z = 0.5. The tendencies are the same. For zero momentum (Fig. |3|) 
there is no distinction between particles and holes, and the two ^-function peaks correspond 
to positive and negative energy solutions. As the S-field is increased, there is a broadening 
in the width and the positions of the peaks are shifted towards the zero temperature value, 
which at p z = is a peak at E = 0. At non-zero momentum (p z = 0.5, Fig. |j) the two 
peaks at B = correspond to a particle solution at E ~ 1.1 M. e and an anti-hole solution 
at E ~ — M. e . In addition there is a continuous part in the interval — p z < E < p z . Also 
in this case the peaks get broader as B increases and eventually there is only one wide peak 
around p z = 0.5 for very strong fields. In vacuum there is an imaginary part of the self-energy, 
describing the decay to a lower energy level due to synchrotron radiation, only for the higher 
Landau levels, but at finite temperature even the particle in the lowest Landau level can 
scatter with the surrounding plasma and this is what causes the imaginary part. 

It should be noticed here that, as already discussed at the end of section 0, the HTL 
approximation is only valid for momenta smaller than the temperature. For larger momenta, 
or stronger fields, the tree level contribution dominates the real part in the Dirac equation, 
but since there is no imaginary part at tree-level it comes entirely from the HTL term, 
which is then not reliable (though small). In the very strong field limit (qB T 2 ), only 
the lowest Landau level is occupied and the imaginary part comes from annihilations with 
the antiparticles |7[|. This term is linear in temperature and does not appear in the HTL 
approximation. 



5 The chiral anomaly 

The classical action for massless fermions is invariant under chiral transformations, but the 
corresponding chiral current is not conserved on the quantum level due to the chiral anomaly 
||. The divergence of the chiral current in 3+1 dimensions is given by 

djg = w*75* = w^^ pCT • (5-1) 
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Figure 4: Spectral function for various magnetic field strengths at the momentum 
p z = 0.5 in the lowest Landau level. For very weak fields there are two 5-functions at 
the positions of the particle and anti-hole poles (indicated by the vertical solid lines), 
and a continuous part in the interval [—0.5, 0.5] which is below the light cone. When the 
field increases the anti-hole peak around E ~ —1.0 disappears faster than the particle 
peak at E ~ 1.3. At qB = 5-OA^e there is only a very wide peak left of the particle 
around E ~ 1.0, which for increasing field is slowly shifted towards E = 0.5 while also 
getting narrower. 



Finite temperature effects do not break chirality and, as a classical action, the HTL effective 
action is still chirally invariant. Since the anomalous term in Eq. Q5.1J) originates from the 
UV-divergent part of the propagator, it is not expected to change at finite temperature and 
Eq. ( |5.1|) is expected to hold. This has also been verified explicitly by several authors [10. 1 1 



In vacuum there is a clear physical picture, related to the IR properties of the fields, of how 
chirality is created by moving particles from the Dirac sea up to positive energy by switching 
on an external electric field [|T2|, |T3[. This picture works because the massless dispersion 
relation crosses the Dirac surface and particle-antiparticle pairs can be created continuously. 
By adding a chirality-breaking mass term the dispersion relation no longer crosses the Dirac 



15 



surface and, in fact, no chirality is produced if the external gauge field varies adiabatically. 
New particle pairs are again created |13| when the variation of the gauge field is rapid with 
respect to the mass of the fermion. 

At very high temperature, it is well known that even massless chiral fermions pick up an 
effective thermal mass, in the sense that the energy of the propagating modes does not go 
to zero for vanishing momentum. From this IR picture it is not obvious how the anomaly 
equation (Eq. ( |5.1[) ) can be fulfilled at finite temperature. What happens, as I show be- 
low, is that as a quasi-particle of one chirality moves along the dispersion relation from a 
particle-like excitation to a hole-like one, it looses its spectral weight, while a quasi-particle of 
opposite chirality gains the same amount of spectral weight. In this way the spectral weight 
is shifted between chiralities (or between states above and below the Dirac surface) without 
any dispersion curve actually crossing the surface. 

Most discussions of the anomaly at zero temperature are performed using a language of 
single-particle excitations. Even though there are stable quasi-particles in the leading HTL 
effective action without any external field, this is far from being the whole picture. There is 
some danger in treating the branches of the dispersion relation as ordinary particles. Each 
of the branches does not have a full spectral weight, and not even the sum of the spectral 
weight for the particle and the hole adds up to 1. With non-zero external field there are 
also imaginary parts, which cannot be accounted for within a quasi-particle picture. We 
shall therefore carry out the calculations entirely in terms of Green's functions and without 
reference to single-particle states. The HTL effective action is, after all, only a way of writing 
a set of Green's functions. 



5.1 The 1 + 1 D anomaly 

In order to see what are the essential parts of the thermal anomaly equation, we shall first 
briefly repeat the standard calculations of the anomaly in 1+1 dimensions using both an 
operator and a Green's functions language. We shall follow ref. |13| very closely in the 
operator formalism. The chiral anomaly at finite temperature for a non-interacting theory 
has been studied by several authors |TUJ, in particular in 1+1 dimensions. While Smilga [TT[ 



gave a physical interpretation of the thermal effect in terms of scattering with particles in 
the heat bath, I take a slightly different approach to reach similar conclusions. Once we have 
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understood the mechanism in 1+1 dimensions it should be easier to see that the same holds 
true in 3+1 dimensions. 

There is an intrinsic consistency problem in treating a time-dependent gauge field, i.e. an 
electric field, and an equilibrium ensemble at the same time. Starting from an equilibrium 
ensemble it will not remain in equilibrium if we switch on an electric field, unless we consider 
an adiabatic limit where the ensemble has time to readjust itself to equilibrium much faster 
than the variation of the gauge field. Even though Eq. ( |5.1| ) is true at the operator level, and 
thus true for any expectation value of the equation, it does not determine the time-integrated 
form (it'tr [p(t')\l/(x)7o75\l/(x)] for an explicitly time-dependent density matrix p(t). The 
fields in Eq. ( |5.1| ) are in the Heisenberg picture and if the density matrix has no explicit time 
dependence it enters only as an initial condition. In an adiabatic limit one could effectively 
take into account interactions with the heat bath by using a density matrix, which at all 
times corresponds to thermal equilibrium. This would then be an explicitly time- dependent 
density matrix and the time-integrated anomaly equation may not fulfil the standard anomaly 
equation. Such a replacement with an effective density matrix is to some extent arbitrary and 
depends on which physical situation is imagined. It is possible to consider that only energy is 
equilibrated by scattering processes, and that chirality is conserved in each process, and thus 
not equilibrated. Or, one can consider the system to be in contact with a heat reservoir with 
which it can also exchange chirality. 

Another possible situation is when we neglect interactions between particles altogether 
and follow the exact time evolution of the non-interacting plasma, after its initial condition 
is given. This is the situation we shall consider in detail. 

The chiral charge has to be defined using a gauge-invariant point splitting regularization 
in the spatial z-direction: 



J Jn, 



(5.2) 



where 

r e 2 
exphf^j 



W = ■ (5.3) 



Including a chirality-breaking Dirac mass the anomaly equation in 1+1 dimensions, integrated 
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over space and over time from to r, reads 

hm(Ql(r)) = ^- f dt [dx e^F^ + 2im f dt{Ql{t)) , (5.4) 
7^0 In Jo J Jo 

where (Ql) is defined as {QD in Eq. (|5~^ ) but with 7075 replaced by 75. With a field operator 

*(*,€■)= /[dpJe*" € '[ti p .(06p.+«- fl ,(*Aj , (5-5) 



where the notation is taken from [ 13 1 , the chiral charge in the massless limit can be computed 

as 

L 



(QUr)) = ~J dp z e-^-^)?(l - (blbj - (dl Pz d. Pz ) 



0(Pz) - o(- Pz 



(5.6) 



The initial expectation values of the number of particles (b p b Pz ) and antiparticles (d}_ Pz d_ Pz ) 
depend on which physical situation we consider, but in any case they should go rapidly to zero 
for large \p z \. In thermal equilibrium, with zero chemical potential, we would for instance have 
{b^ p b Pz ) = {d}_ Pz d- Pz ) = (exp[/3\p z \] + l) -1 . It is thus only in the vacuum part that the point 
splitting is needed. It follows that only the vacuum part can depend on A z (t) and the chirality 
production is, therefore, independent of the initial thermal condition. If, on the other hand, 
we consider a situation where the particles relax rapidly to thermal equilibrium, so that the 
distribution of particles with quantum number p z is determined by the energy of the states 
after switching on the A z -field, then we should rather use {b p b Pz ) = (exjp[/3E Pz ] + l)" 1 , where 
E Pz = \p z \ — siga(p z )eA z (t). In this case also the thermal part of Eq. ( |5.6|) depends on A z (t), 
not directly through the anomaly, but from the interaction with the heat bath that we assume 
to be present in order to maintain thermal equilibrium. 

Since in the (3+l)-dimensional case we want to avoid the use of single particle states 
as in Eq. ( p.5| ), we shall now see how the above calculation can be performed using Green's 
functions. The field expectation values can be related to the time-ordered Feynman Green's 
function via 

(* a (OMri)) = -iS F (r l ,t) ba \ io>Vo . (5.7) 

In a time-dependent background field, where energy is not conserved, it is not obvious what 
the correct e-prescription for a time-ordered Green's function should be. In this simple non- 
interacting case we can, however, compute everything explicitly and, starting with the vacuum 
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part, we find 



T*(£)*(77)>° = J[d 2 p}ex V -ipoto + Wztz + *15 f° dtA z (t) 



exp 



WoVo ~ iPzVz + il5 / dtA z (t) 



</<> 



(5.f 



where the usual Feynman prescription p — > (1 + ie)po should be used in j. All complications 
from the external A z field is thus taken into account in the phases of the diagonalizing wave 
functions, and the Feynman prescription is unchanged. The significance of po depends on the 
representation of the wave functions we use to diagonalize the propagator, and here po is the 
initial energy of a particle in a state labelled by p z before the A z field is turned on. The actual 
energy of that state then varies as po — xA z (t), where \ is the chirality of the state. 
With this form of the propagator we obtain the chiral charge from 



(*(07o75* fo)>° = [d 2 P ]e 



ip{£-rj) 



-i f io eA z (t)dt 

e j to 



i f ( ° eA z (t)dt 

e Jr >o 



Po~Pz + itPo Po+Pz + itPo 



(5.9) 



Zo>Vo 



The condition £ > tells us that the po contour must be closed in the upper half-plane. 
The poles in Eq. ( |5.9|) give two ^-functions in p z in the standard manner. The rest of the 



calculation can be found in |13j and the result is 



(Q1(t))° 



L 
2^ 



dp z e 



r(p*-eA«) s 



[0(Pz) - 0(- Pz )} 



7^0 



2tt 



-2eAJr) 



(5.10) 



in accordance with Eq. (|5.6|) . From this exercise we learn that in the massless case, where 
we can find an explicit basis diagonalizing the propagator, the mathematical mechanism that 
gives us the correct anomaly is that the poles in p cross the real axes when p z = 0. 

In the equilibrium real-time finite temperature formalism the free propagator can be writ- 
ten as in Eq. (|4.2|) , but since we compute a one-point function we only need the 11-part: 



iS&p) = *S° F (P) - /f(po) (iS° F (p) - *S° F *(p) 



(5.11) 



where /f(po) is the thermal distribution function. The problem we have at hand is not one 
of an equilibrium, but as we saw in Eq. ( p.8[ ) the time dependence can be entirely absorbed 
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in the phases of the wave functions. Since in this basis po has the meaning of the energy of 
the initial state, the thermal version of Eq. ( |5.9|) is obtained by the substitution 

2irsign(p )f F (p Q )A x (p ,p z ) = 2irsign(p Q )f F (p )5(p -xPz) , (5.12) 



Po - XVz + itPo 



where /f(po) is the initial particle distribution. With this propagator the thermal contribution 
to the anomaly is given by 
L 



(Ql{r)Y = ^J dp z e-^- eA ^ 2 J™ dp (f F (p ) + f F (-p Q ] 



A H (po,p z )-A L (p ,p z ) 



(5.13) 



which agrees with the thermal part of Eq. ( |5.6| ). Since the spectral functions are rapidly con- 
vergent in p z for fixed po the point splitting 7 can be sent to zero before doing the integrations. 
There is no A z (t) dependence left, which shows again that there is no thermal correction to 
the anomaly. 

5.2 The anomaly in 3+1 dimensions at high temperature 

The HTL effective action is chirally invariant even though there is a mass gap in the dispersion 
relation. We, therefore, expect that the chirality produced in vacuum cannot be undone by 
chirality-conserving interaction with the thermal heat bath. From a mathematical point of 
view we can argue that since the anomaly equation is true at the operator level it must 
remain true in whatever average we take, including a thermal average. We shall see that 
this is correct by explicitly calculating the chirality production in 3+1 dimensions within the 
HTL approximation. The only thing we need is an explicit expression for the propagator. The 
fermionic part of the HTL effective action is simply related to the inverse of the propagator 
by 

C i HTL =^(x)S- 1 (x,y)^(y) . (5.14) 

In order to write down the propagator itself we have fixed the boundary conditions when in- 
verting the kernel of Eq. ( |5.14j ) . At finite temperature the standard inversion gives Eq. ( |5. L 1[ ) , 
but with 



iS° F (x,y) = (T[*(x)*{y)]) = {x 



u-ri 



y) • (5.15) 
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The e-prescription for time ordering can be obtained by comparing the present calculation 



with the explicit calculation in section |5.1| . We shall here use similar, explictly time-dependent, 
wave functions to diagonalize the propagator. 

The (3+l)-dimensional anomaly equation at finite temperature, in a classical background 
field, is given by 



e 2 



(d^(x)Yl 5 V(x)) = j^e^F^ . (5.16) 

We shall compute the left-hand side of Eq. Q5.16j ) from the propagator in Eq. ( |5.15| ) in a back- 
ground consisting of orthogonal E and B fields. Scattering with external thermal particles, 
described by the thermal part of Eq. ( |5.11| ), is discussed at the end. In the massless vacuum 
case we saw that the essential mechanism of generating the correct anomaly was that the 
poles of the propagator crossed the real axis when p z changed sign. This does not happen at 
finite temperature due to thermal masses. 

First, we need to diagonalize the HTL Dirac equation in the presence of a S-field in the 
z-direction and a parallel E-fie\d. We choose the gauge = (0, 0, 0, y4 3 (£ )) for the electric 
part of the background field. The diagonalization can be done in the same way as in section ^ 
but instead of Eq. ( |2.8| ) we use the basis 



(£K) = exp 



V z fso I 

-ipo^o + ie— J A z (t)dt + ip z £ z + ip y £ y + i \p x £ x + 



(5.17) 



but the eigenvalue of v • II remains v ■ p. In this way the propagator can be calculated 
exactly, even for non-adiabatic electric background fields, just as in the massless case at zero 
temperature. The prescription for the time-ordered Green's function is again po — > (1 + ie)po, 
in this basis. We write the relevant trace of the propagator as 

4 

tr S F (y,x)>y >y 5 = J2 J2 (v> a\K,i)[S F (n)^^} i:j {K, j\x, a) , (5.18) 

K i,j,a=l 

and use the basis 



■ 



(5.19) 
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where is a set of 4-spinor base vectors, which can be taken to be u£ = 5^,. For n = there 
are only two states, vS 1 ' and u^ 3 \ the others being identically zero. With the Dirac operator 
in Eq. ( |2.16| ), diagonalized in the spatial quantum numbers, we obtain for the higher Landau 
levels (n > 0): 

tr5 F ( K ) 7 o75 = tr J D^ 1 (/ t )-tr J D^ 1 (K) . (5.20) 

In the lowest Landau level the matrices Drl are not invertible, since there are only two states 
in total. Its explicit form then is (k = {n = 0,p y ,p z }): 

J-oo p — E + tep 



with A^ L (E,p z ) given by 

1 



Po-p z - M 2 e ((u ) K0 - (u z ) Ko ) 



°° dE ^MEl^L (5.22) 
-oo po - E + lepo 



and Alll(E,p z ) = A^ LL (E, —p z ). We shall start by computing the contribution to (Q5) from 
the lowest Landau level. This is, in fact, the only part that contributes, as we shall see later. 
Taking the point splitting only in the z-direction we find 

(Q 5 >lll = -iV J deA(&) / [^o]fe][^](0|K )tr5 F («:o)7o75(^o|Oe ieAz6 ■ (5.23) 

Using 

r eB 

J[d Py ]I n ^ x ) = — , (5.24) 

and the fact that the po-contour should be closed in the upper half-plane, the produced chiral 
charge reduces to 

(Q 5 >lll = -^f / dp z e-i^-^f [ W {p z ) - W{-p z )} . (5.25) 

This equation has a clear resemblance with Eq. ( |5.10| ). The function W(p z ) is defined by 

W(p z )= / dEA* LL (E,p z ) . (5.26) 
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2 F 



Po 




Pz 

Figure 5: A density plot of the spectral function A^ LL (po,p z ) for qB = M. 2 e . Even 
though there is no crossing of the Dirac surface, the spectral density goes continuously 
between positive and negative energy states when p z decreases. 

It is the spectral weight for the right-handed positive energy solution in the lowest Landau 
level. For very large \p z \ there are no collective excitations, such as holes, but only the 
standard particle solution. With our convention that the 5-field points in the positive z- 
direction, we find that for p z > 0, A^ LL (E,p z ) is concentrated at a positive energy particle 
solution at E = \p z \, while for p z < it is peaked on a negative energy antiparticle solution 
at E — — \p x \ (see Fig. ||). Thus, we have 

lim W(p z ) = , lim W(p z ) = 1 , (5.27) 

and all derivatives of W(p z ) vanish for large \p z \. It can then be shown that 

VeB r r l e 2 

(Q 5 (*)>lll = -^~tzA z = J d 3 x J dt'—e^F^ . (5.28) 
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We note that this agrees with the anomaly in Eq. ( |5.16| ) for the particular background field 



configuration that we have chosen. When it comes to the higher Landau levels it turns out 
that the two terms in Eq. ( |5.20 ) are separately well convergent for large p z , as opposed to 



Eq. ( |5.2ip , which needs the point splitting in order to be well defined. We can, therefore, 
change p z — > — p z in Dl(k), after which the sum cancels when 7 — > 0. 

Scattering with particles in the thermal heat bath is taken into account in the same way 
as in 1+1 dimensions. The thermal part is UV-convergent and, exactly as in Eq. (|5.13|) , it 
has no time dependence. 

6 Conclusions 

The main computational part of this paper is the diagonalization of the fermionic part of 
the Hard Thermal Loop effective action in the presence of a constant background magnetic 
field. This makes it possible to write down the explicit expression for the spectral function 
of fermions and to see how it depends on the magnetic field strength. We find that, starting 
from weak fields, the spectral weight moves from the standard particle and hole solutions at 
high temperature over to the vacuum Landau levels for very strong fields. 

It has been recognized in the literature that it is difficult to reconcile the standard picture 
of level crossing as a mechanism for anomalies, with the thermal masses of fermions at high 



temperature ||11|| . Using the exact spectral function, I have shown in this paper how the 
spectral weight can move continuously between chiralities, when a background electric field 
is switched on, without the dispersion relation ever crossing the Dirac surface. The differ- 
ence from the vacuum is that the spectral weight on a dispersion curve varies continuously 
between zero and one in the HTL approximation, while it is always exactly one in vacuum. 
The conclusion is that the anomaly equation remains valid at high temperature, even after 
taking interactions into account. It is however possible to obtain different production rates by 
coupling the system to an external heat reservoir, which means effectively using an explicitly 
time-dependent density matrix. The chirality production in that case would then depend on 
the exact experimental setup, and I have not discussed this possibility in any detail. 

Even though the main problems formulated in this paper have also been solved here, there 
are some related issues that still call for solutions. One problem is to extend the analysis 
of the dispersion relation to non-Abelian gauge bosons, but this turns out to be far more 
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complicated due to self-interaction. Another problem of a certain interest is to see how this 
anomaly mechanism fits into the language of index theorems, which has shown to be useful 
for the anomalies at zero temperature. 
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A Normalization and causality 

The thermal expectation value of the canonical anticommutation relation for the fermionic 
fields 

({*(*), *t(y)})=C(x-y) (A.1) 

should vanish for space-like x — y and should be equal to a ^-function in x — y when xq = xjq. 
These are basic requirements of the fundamental fields which we do not expect to be violated 
by the HTL approximation. In terms of the spectral function the normalization condition, 
derived from the equal-time commutator, becomes 

[°° dEA(E lP ) = ^- f°° dE(S(E-te,p)-S(E + te,p))= l0 . (A.2) 



There are two ways of showing the validity of Eq. ( |A.2j ) for the propagator in Eq. ( |4.5| ), 
analytically and numerically. First we present the analytic proof. Consider the contribution 
from the advanced and the retarded propagators separately. For a given p they are analytic in 
the lower and upper half-plane, respectively ||14|| , and the integration contour can be deformed 
to suitable arcs at infinity. The integrals along those arcs do not vanish, but the HTL 
corrections go to zero. The expression for a free fermion can then be used and it trivially 
satisfies Eq. ( |A.2[ ). In this way it is clear that the only important ingredients for the spectral 
sum rule to be fulfilled is the analyticity away from the real i?-axis, and that the correction 
goes away for large complex E. It is only the 70-part of the propagator in Eq. (|4.5| ) that 
contributes to Eq. (|A.2|) since the other parts of Eq. (|4.5| ) decay too fast on the arcs at infinity, 
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and they are also antisymmetric in E. It is, therefore, common to write the normalization 
condition as 

/ dE-tr[ 7o A(E,p)} = l . (A.3) 

J — oo 4 



The other way to check Eq. (|A.3| ) is by a direct numerical calculation. The residues at 
the poles above the light cone (E > \p\) have been computed by several authors |5], and it 
is well known that they do not add up to 1 for \p\ > 0. It is straightforward to calculate the 
integral in Eq. (|A.3| ) below the light cone, and it turns out to make up for the missing part, 
as expected. 

The causality condition, i.e. that Cp{x) vanishes for space-like x, has been discussed in 



the HTL approximation in [16], where 

dEd 3 p 



C{x) 



(2tt)* 



exp[-i(Et - px)]A(E, p) 



(A.4) 



was calculated numerically for gauge bosons. I shall here give an analytic demonstration that 
the commutator indeed vanishes outside the light cone, starting with the part of Eq. (|A.4|) 
proportional to 7 , C (x) = tr jj C(x). Assuming that tr [ r y A(E,p)] only depends onp = \p\, 
we can perform the angular integral and obtain 

ip\x\ _ 

11- 

L 4 



C (x) 



(27r) 2 J-oo 



oo poo 

dE / dpp 2 e 
J o 



2-iEV 



vp\x\ 



-ip\x\ | 

tr l- l0 A(E,p)} 



1 d r°° , 1 

dEdpe-* Et+ ^tri- l0 A(E,p)} . 
4ii 2 \x\ d\x\ J-oo 4 



(A.5) 



In order to more easily see the analytic structure, we change the variables to u (= E +p) and 
v (= E — p): 



C (x) 



1 



d r°° 



8n 2 \x\ d\x\ 



du dv exp 



\<t-\x\)- l -v(t+\x\) 



1 

x tr 

2ni 



-70 yS(u —ie,v — ie) — S(u + ie, v + ie] 



Let us first see how it works for a free massless scalar, where 

1 1 



S(E-ie,p) 



E 2 —p 2 —iEe u—~€V — \e 



(A.6) 



(A.7) 
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The integrals over u and v factorize and we can use 



du = 2ni9(a) , (A. 

oo u — ie 



to show thatf] 



c (x) = -^^L[e(-t + \x\)e(-t-\x\)-e(t-\x\)e(t+\x\)) 

-[6(t-\x\)-5(t+\x\)] = -^-sign(t)5(t 2 ~\x\ 2 ) . (A.9) 



4rc\x\ 2ty 
From this follows also the canonical commutation relation for scalar fields: 

= dvo Co(x - = i5^\x - y) . (A.10) 

The basic reason why the commutator vanishes outside the light cone is the occurrence of 
^-functions coming from Eq. ([A. 8 ) and this follows from the property of analyticity in the 



correct half-plane. Since Sq(u — ie,v — ie) is analytic in the lower half-plane for both u and 
v, and it vanishes fast enough for large arguments, the integration contours can be closed in 
the lower half-plane for positive t — \x\ or t + \x\, respectively. We conclude that for any 
S(u — ie,v — ie) with the correct analyticity properties we have, using 9(a)9(/3) = 9(aj3)9(j3), 



du dv exp 



± u (t-\ x \)- % -v{t + \x\) 



S(u-ie,v-ie) =9(t 2 -\x\ 2 )F(t,\x\) , (A.ll) 



for some function F(t, \x\), at least for non-zero t 2 - \x\ 2 . There can be other singularities 
right on the light cone. A similar argument applies to Sq(u + ie, v + ie), leading to another 
factor 9(t 2 - \x\ 2 ). 

Let us return to the HTL propagator and see if it has the correct analyticity properties. 
It is given by Eq. (|4.5|) , and in terms of u and v we have, for E — > E — ie: 

2M 2 e , fu-ie 
In 



u 2 — v 2 v r — k 



AM 2 (. u + v , (u-ie\ , , 



2 Note that there is a sign error in the revised version of ref. ||i 
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The cuts from the logarithm start slightly above the real axis and the branch cuts remain 
in the upper half-plane. There is a potential singularity at u = v, but it can be shown that 
both s(u,v) and r(u,v) have power series expansion around that point. Thus, none of these 
functions have any non-analyticity in the lower half-plane. Then we only have to verify that 
the denominator in Eq. ( |4.5| ) does not have any pole in the lower half-plane. This should be 
done for each of u and v keeping the other one real, and it is not very difficult to do this 
numerically. With a fine enough grid one can demonstrate that there are no singularities 
away from the real axis. When deriving Eq. (|A.5 ), we assumed that the spectral function 
only depended on \p\. The HTL spectral function also has a term proportional to 7p, but 
it can be rewritten as i~y ■ V acting on a rotationally invariant function, so that it does not 
really affect the above reasoning. 
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